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Abstract: We investigate factorisation at small x using a variety of analytical 
and numerical techniques. Previous results on factorisation in collinear models are 
generalised to the case of the full BFKL equation, and illustrated in the exam- 



D . pie of a collinear model which includes higher twist terms. Unlike the simplest 

collinear model, the BFKL equation leads to effective anomalous dimensions con- 



taining higher-twist pieces which grow as a (non-perturbative) power at small x. 
While these pieces dominate the effective splitting function at very small x they 
5^ \ do not lead to a break-down of factorisation insofar as their effect on the predicted 

scaling violations remains strongly suppressed. 
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1. Introduction 

The standard approach to calculating any quantity in Deep Inelastic Scattering (DIS) 
or Drell-Yan type processes is to assume factorisation, based on the operator product 
expansion. This allows us to consider parton distributions which evolve with scale 
as 

?^^,M')UQ^Ho[l,), (1.1) 

where u is the Mellin moment of the parton distribution. Cross sections are then 
expressed in terms of the parton distribution at scale Q^, multiplied by a coefficient 
function, which is itself a function of only uj and asiQ"^), plus corrections of order 

i/g=^". 

A few years ago however, it was suggested in [1] that at small x (small u) 
diffusion might place a limit on factorisation: namely that beyond a certain value 
of X factorisation would break down. The argument was the following: at moderate 
values of x parton evolution in Q^ can be associated with a chain of emissions ordered 
in Q^. Corrections to factorisation can arise from chains which evolve from some 
starting scale Qgtart ^P ^'^ ^ scale Q^, evolve down to a non-perturbative scale A^ and 
then evolve back up to scale Q^. At moderate x, such contributions are suppressed 



by a power of h? /Q"^. At small x however, this process of going up and down in scale 
becomes an essential element of the evolution, and is no longer suppressed by powers 
of Q^: so it would seemed that factorisation might break down at the values of x 
such that this diffusion enters into the non-perturbative region. More recently this 
has been used as an argument for parameterising and fitting the small-x behaviour 
of splitting functions, rather than trying to calculate them from first principles [2] . 

In [3] plausibility arguments were presented suggesting that there would be no 
such breakdown. Physically speaking, perhaps the simplest explanation is to say 
that evolution that, while going down in x drops below one's starting scale Qstart) 
contributes to the small- a; part of the gluon distribution at the starting scale. There- 
fore it should not be included in the anomalous dimension in order to avoid counting 
it twice, so the anomalous dimension contains only evolution which remains above 
the starting scale, ensuring that it is perturbative. 

However this explanation, and the arguments in [3] , are not entirely satisfactory 
because they are not able to estimate the size or functional dependence on x of any 
violations of factorisation. 

In [4], a model was studied which was qualitatively similar to BFKL in that it 
had diffusion, and which correctly reproduced collinear evolution at all orders. It 
was possible to show that this model leads to exact factorisation, namely that there 
were no 1/Q^ corrections at all. Thus diffusion on its own is not incompatible with 
factorisation. 

This collinear model had two particularities: the branching kernel, when con- 
verted to 7 space (the Mellin transform variable conjugate to transverse momentum) 
contained poles only at 7 = and 7 = 1, i.e. the leading collinear and anti-collinear 
poles. The BFKL equation instead has poles at all integer values of 7. The second 
particularity was related to the scale of a^'. in the collinear model the scale was cho- 
sen to be max(A;^, k' ), where k and k' are the exchanged transverse momenta before 
and after the branching. In the BFKL equation, the next-to-leading corrections seem 
to indicate that the correct scale is rather q = \k — k'\ [5] which has the property 
that it can go to zero even when both k and k' are large. 

Thus there is a need for a study of factorisation within the BFKL equation. 
The complications mentioned above make this quite difficult to do analytically. We 
therefore adopt a two-pronged approach. In section 3 we consider an extension of 
the collinear model with additional poles at 7 = — 1 and 7 = 2. Some of its gross 
features can be deduced analytically and can be expected to carry over to the full 
BFKL equation. 

In section 4 we then present a simple numerical method which for the first time 
allows the extraction of the exact (effective — see below) anomalous dimensions in the 
full leading logarithmic BFKL equation, with all-order running coupling corrections, 
allowing a detailed study of factorisation not limited by the simplifying assumptions 
of any given particular model. 



2. Understanding factorisation 

We shall start off by discussing factorisation for the non-integrated gluon distribu- 
tion, or equivalently for the gluon Green's function. Conceptually it it somewhat 
simpler — the generalisation to the integrated gluon distribution is then relatively 
st raight forward . 

We shall discuss factorisation by considering a number of models. We study a 
gluon Green's function G^(t, tg), t = InQ^, in analogy with the unintegrated gluon 
distribution in the proton, and consider as the non-perturbative aspects of the 
problem the value of to (the lower hard scale of the process) and the regular isat ion of 
the running coupling in the infra-red. The latter will be represented by some scale i 
at which the coupling is cut off (we could equally have chosen to freeze the coupling 
at that scale). The small-x properties of the Green's function will depend on both to 
and t. More specifically the position of the leading pole in the u;-plane will depend 
on the value of t and its normalisation on both t and to. 

If factorisation holds then the (non-integrated) effective anomalous dimension 
defined by 

~ /^ 1 dG^{t,to) 

should be independent of both t and to at least to within higher-twist terms sup- 
pressed at least as e~*. 

This is not quite a sufficient condition: indeed if the anomalous dimension con- 
tains higher-twist pieces which grow as a sufficiently large power of x then these 
could dominate the scaling violations making them impossible to predict. Since the 
Pomeron singularity will affect G^^ anyway, yielding a large power behaviour of this 
kind, the problem arises of understanding whether approximate factorisation is still 
preserved in the small-x region. 

Since we are interested in pieces enhanced at small-x we need to understand the 
singularity structure of the anomalous dimension. There are two possible origins 
for singularities of j^. One possible origin is that at some u, dG^/G^ contains a 
non-factorisable singularity, the other is for G^ to be zero while dG^ is finite. 

Now, the simplest collinear model of [4] provided us with a mechanism by which 
the Pomeron singularity in G^j can be consistent with exact factorisation: the leading 
singularity of % was found to come from a t-dependent zero of G^, leaving the 
Pomeron factorised away. We shall see in the following that this basic mechanism is 
still at work in collinear models with higher twist terms [6] and in the BFKL equation 
itself, thus suppressing the non-factorisable singularities. 

2.1 Recalling the structure of the collinear model 

The collinear model of ref. [4] gave a powerful handle for the study of anomalous 



dimensions. It was described by an equation of the form 

ujG^it,to) = 5it - to) + fdt'K2it,t') G^it'M), (2.2) 

whose kernel is 

K^it.t') = a,{t)Q{t - t') + as(t')e(t' - t)e-^''''\ (2.3) 

In 7-Mellin transform space with respect to Q^ = exp(t), the leading order part of 
the kernel gives the following characteristic function: 

X'-'^'^I) = - + -r^ , (2.4) 

7 1-7 

for which reason it is also referred to as a 2-pole model. It has the convenient 
property that it can be expressed in terms of a second order differential equation, 
whose solutions have the following factorised form: 

G^{t,to)=F^{to)F^{t), t>to, (2.5) 

where F^ and F^ are the linearly independent solutions of the homogeneous equation 
(eq. (2.2) without the (5-function term), that are regular to the left (negative to) ^-^d 
to the right (positive t) respectively. 

The fact that F^{t) is regular at large t means that it is independent of t for t > i. 
It is free of singularities in cu-space, as illustrated in figure 1. All non-perturbative 
dependence, in particular the poles governing the high-energy behaviour of G, is 
contained in F^, also illustrated in figure 1. 

One sees from these plots that zeroes of G^ can arise from both F^ and F^. ^ 
However since the t-dependence lies entirely in F^, the to-dependent, non-factorisable 
zeroes of F^ lead to zeroes of both G^^ and of its t-derivative. Therefore they do not 
lead to divergences of the anomalous dimension. Only the t-dependent zeroes of F^ 
lead to poles in 7(0;), dtF^{t) not usually being zero when F^{t) is zero. 

This fact has the consequence mentioned before: the Pomeron singularity, present 
in F^{to), does not occur in the anomalous dimension, which is singular at ujc{t), the 
leading zero of F^ (while we refer to the leading zero of F^{to) as cjo(to))- 

3. The 4-pole collinear model 

One of the main differences between the BFKL equation and the collinear model just- 
discussed lies in the presence in the BFKL equation of poles of xil) ^^ ^-H integer 

^Note however that F^ does not always have zeroes: their presence depends on a variety of factors 
such as the relative sign of successive divergences, which are determined by the non-perturbative 
parameters. Only in the situation of to 3> < is at least one zero guaranteed, just to the left of the 
leading divergence. 
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Figure 1: Illustration of the kind of the w-dependence that arises in the right- and left- 
regular solutions of models based on second-order differential equations such as the two-pole 
collinear model. (The small-u; parts of the curves have been left out in order to improve 
the overall clarity of the plots) . 



values of 7. The poles beyond 7 = and 7 = 1 can give rise to higher-twist effects, 
in which we are particularly interested in this article. 

A useful stepping stone to the full BFKL equation is an equation containing the 
first subleading higher- twist parts, namely poles at 7 = — 1 and 7 = 2. The equation 
for the Green's function is 



ujG^it, to) = Sit - to) + / dt'K^it, t') G^(t\ to). 



(3.1) 



where the kernel is 

K^{t,t') = as{t)Q{t - f) {l + e-(*-*')) + a,{t')Q{t' - t) [e-^''-'^ + e-2(*'-*)) 

-^a.{t)5{t-t'). (3.2) 

The (5-function term is included so that the value of Xo(l/2) be the same as for the 
2-pole coUinear model: 

x'-''°^^(7) = T^ + - + T^ + ^-t (3-3) 

1 + 7 7 1 — 7 2 — 7 6 

The physics content of this formulation is quite similar to the leading and first- 
subleading collinear and anti-collinear parts of the BFKL equation. A slight differ- 
ence exists in that the higher-order higher-twist pieces of the BFKL equation also 
contain a non-local dependence on as- Their inclusion would considerably increase 
the difficulty of the analytic treatment. 

This model as it stands can be expressed as a fourth-order differential equation 
(or two coupled second order differential equations, etc.), using the same approach 
as in [4], as discussed in detail in the Appendix. It is thus similar to the diffusion 
models with higher twist terms [6]. Without having to write down the full equations 
in detail, we describe here some important properties of its solutions. 

There are four linearly-independent solutions of the homogeneous equation. We 
denote them by F^'^, F^'^, the leading and sub- leading twist right-regular solutions, 
which for large t go as a constant and e~* respectively; and by F^'", F^'^ for the two 
left-regular solutions. For t > t^ the Green's function has the expression 

G.(t,to) = Fi'"(to)Fi^'°(t) + Ft'\u)F^^\t) . (3.4) 

as expected from [6], and proved in the Appendix. There we also show how to de- 
termine the F-^'s so as to satisfy the boundary conditions at t = to- The solution for 
t < to can then also be obtained, by a simple exchange of variables. As in the 2-pole 
model, the right-regular solutions are independent of the non-perturbative parame- 
ters of the problem. The Pomeron singularity and the non-perturbative dependence 
all enter into the left-regular solutions. 

It is convenient to divide the solution into its leading and higher-twist parts, 
Gj^ and G^, corresponding to the first and second terms of (3.4) respectively. For 
large t ^ to,t, we have that G^ is strongly suppressed (by a relative amount ~ 
e~*) compared to G°; G° should be qualitatively similar to the solution of the 2- 
pole model, having zeroes of perturbative and non-perturbative origin. We know 
therefore, from the start, that the violation of factorisation is uniformly of higher 
twist, just because of the additive nature of (3.4). 

On the other hand, one can look at the same problem from the standpoint of 
the singularities of the anomalous dimension. Let us concentrate on the case with a 



discrete spectrum (running coupling with cutoff), in which we have to look for zeroes 
of Gi^. The fact that G^ is strongly suppressed means that the distribution of zeroes 
of G^^ is determined essentially by the positions of the zeroes of G° . But in general 
where G° is zero, G]^ will be non-zero, causing the zeroes of G^^ to be slightly shifted 
compared to the zeroes of G° . The size of this shift will be related to the relative 
sizes of G^, and G,",, i.e. it will be of the order of of e^*. 

The effect of the shift will be different according to whether the zero is of per- 
turbative or non-perturbative origin. In the case of the leading zero of perturbative 
origin, then the leading perturbative pole of the anomalous dimension is shifted, and 
its normalisation changes — both the effects are of order e~* (the relation between 
relative changes to the normalisation and the position of the divergence depends on 
the nature of the t and cj-dependence of G^, and so is difficult to predict). 

The consequences of the rightmost (NP) zero u;o(to) being shifted are some- 
what more interesting, and possibly dangerous for factorisation. We recall that 
in the two-pole model, the non-perturbative zeroes did not lead to divergences 
of the anomalous dimension because the zero was in F^ and so the derivative 
dtG{t,tQ) = F^{tQ)dtF^{t) was also zero when G was zero. This is no longer true 
when higher twist terms are present, and a singularity around ujq is expected. The 
exact value of uq, and indeed even the existence of this zero, depend rather subtly 
on the values of i and to- However, for to significantly larger than t one expects that 
the zero exists, and is driven by the Pomeron term in F^{to), so that cup — c^o ~ e"*" 
is rather small and ujq may be leading. 

To see in detail what happens in the 4-pole model if uq is leading, we approximate 
G°(t,to) around its zero by ^0(6^ - uJo)F^f{t), where J^q = (9^F^'"(to)U=a;o, and 
G^(t,to) by 3\riF^'^(t). It is useful to define partial anomalous dimensions separately 
for F^fit) and F^i^t): 



F^^\t) oj 

where we have neglected terms of order [a^juS)^ in the expansion of the anomalous 
dimensions. Thus F^'^(t) is suppressed compared to F^f(t) by an amount 



^u.W = ^Ru.. = -1 + — + ••• ' (3-6) 



— ^ — oc exp 
F^fit) 



dt' ilUt) - ^^t)) 



-t + 0(a.{t)/a;2) ^ /g^-j^x 



The rightmost zero of G^(t,to), uJq, is therefore shifted compared to that of F^'^(t) 
by an amount 

u'.-uoc:^ - ^'^^y ^ e-* + °(-/-o) . (3.8) 



At Uq the t derivative of G^^ is non-zero, leading to a pole of the anomalous dimension, 
whose residue R{uj'q) is 

RH) = (7io(^) - 7°o(^)) J^l^'o, = (1 + KV^')) K - ^o) (3.9) 

So we have that the residue of the leading NP pole is the same as the shift of the 
zero and both are higher twist. In the case of the BFKL equation (section 4) one can 
test to see if this remains true, in order to verify that the same mechanisms are at 
work there as in the 4-pole model. ^ 

Since uq, which may approach up, is likely to be to the right of the perturbative 
pole Uc, at small x this higher-twist non-perturbative contribution will dominate the 
anomalous dimension. At first sight this might seem to have worrying implications 
for the prediction of small-x scaling violations. But in the present model we know 
that this is not the case, because of the validity of (3.4). From the point of view of 
t-evolution we are saved by the fact that at small x, in the convolution of the splitting 
function with the non-perturbative input distribution (which grows as ujp > ujq), the 
higher-twist part of the effective splitting function gives a contribution of order 



UJp — UJq 

while the perturbative contribution is of order 

Up 

both growing as x""^^. We thus recover in the latter contribution the Pomeron part 
of the leading twist term in (3.4), that we know to be factorised. Thus there is always 
a value of Q^ such that the higher-twist corrections can be ignored at all x. 

A small point worth bearing in mind is that our analysis so far has always been 
for the anomalous dimensions related to unintegrated gluon distributions. In practice 
one is more interested in the anomalous dimensions of the integrated distributions. It 
turns out that their properties are very similar: this is because the n-pole models can 
be expressed in terms of n coupled linear differential equations, and the unintegrated 
and integrated gluon distributions are simply different linear combinations of the 
components of the equations. 



^We note that for the BFKL equation including next-to-leading corrections, the relation between 
- loq and the residue of the pole can be modified by pieces of relative order a^. 



4. The BFKL equation 

We shall study the leading-order BFKL equation including a running coupling, 

g(x, k, Qo) = Si^/Ql - 1) + /" ^ / ^a,{q') 

Jx ^ J '^Q 



X 



A;2 

g(x, l/c - gl, Qo) - 0(/i; - g)S(a:, /c, Qo) 



|A;-g1 



(4.1) 



For notational convenience we have switched to using transverse momenta A;, Qo? 
rather than the logs of their squares t,to. The choice of the emitted transverse 
momentum as the scale for the running coupling, as{q^)., is suggested by the form 
of the NLO corrections to the kernel [5]. A normal DGLAP gluon distribution is 
expressed in terms of S(a^, fc, Qo) via /c-factorisation: 

xgix,Q') = J ^S{x,k,Qo). (4.2) 

In practice rather than solving the integral equation (4.1), it is easier to solve the 
related differential equation 



d9ix,k,Qo) f d'q 
(9 In 1/x 



f d^q 

/ ^as(g^) 



k"^ 

-S{x, \k- q\,Qo) - <d{k - q)S{x, k, Qo) 



\k — g|2 

(4.3) 

with initial condition S(l, k, Qo) = S{k'^/Ql — 1). 

4.1 The extraction of anomalous dimensions 

Naively to obtain the effective splitting function, one would determine S and then 
solve for the function Pgg,cs{z) such that 

xdinQ2gix,Q^) =x —Pgg^esiz)g{x/z,Q'^). (4.4) 

Jx ^ 

However such an approach turns out to be subject to considerable numerical insta- 
bilities. The reason is that any method of solution for S introduces small errors 
(typically of the relative order of 10^^-10^'^). When carrying out the deconvolution 
it generally turns out that P{z) for small z contributes only a small amount to the 
scaling violations {xg{x,Q^) grows as x~'^ whereas perturbatively, P{z) grows as 
x~'^'= and uJc ^ ujf). When x is such that the small error on S{x) is of the same order 
as the contribution to the scaling violations from P{x), then we no longer have a 
handle on the splitting function. 

A solution is to choose an inhomogeneous term such that g{x, Q'^) is independent 
of X. Then, for a given x, the convolution (4.4) is dominated by small z's and small 



errors on g(x) are no longer amplified when translated to P{x). We introduce ?'(x, k) 
as being the unintegrated gluon distribution which, integrated, gives xg{x, Q"^) = 1. 
It satisfies the equation 



d9{x, k) 
Slnl/x 



f{x,k)6{e/Ql 



^/?^^<'^' 








k^ 

— ^{x, k- 

_k-q\^ ^ 


-q\)- 


-Q{k- 


- q)3'{x, k) 



(4.5) 



where 5'(1, k) = 6{k'^/QQ — 1) and f{x) is given implicitly by 



/(^) 



d^A; f d^g 



n 



A;2 



nq^ 



Ois{q^ 



e 



\k-q\ 



-J{x,\k - q\) - e{k - q)'J{x,k) 



(4.6) 



It is trivial to verify that this leads to xg{x, Q"^) = 1. The form of f{x) depends on 
the Q^ value at which we intend to consider the splitting function and on the initial 
scale Qq. Using (4.4) it is now simple to obtain the effective splitting function: 



^_ d3^{x,Q) 



(4.7) 



This turns out to be numerically stable, at least until xPgg^eS becomes comparable 
to the inverse of the machine precision. 

Equation (4.5) is solved by discretising 3^{x, k) uniformally in In k space, and 
then applying standard Runge-Kutta techniques for solving the resulting matrix 
differential equation. This method has the advantage over other potentially faster 
methods, such as a representation with a basis of Chebyshev polynomials [7, 8] , that 
it quite easily accommodates the large variations that arise in the value of 3^{x, k) 
(and errors are at worst of the relative order of the discretisation interval). 

For cSs we take the asymptotic formula, i.e. 



a^{q^ 



Q{q' - Q') 
61ngVA2 ' 



(4. 



where h = 11/12 (we work with zero fiavours). The cutoff at small momenta corre- 
sponds to the regularisation prescription used in the previous section for the collinear 
model. 

4.2 Results 

This section has two aims. Firstly to demonstrate that for the BFKL equation the 
splitting function is a truly perturbative quantity, and that any non-perturbative 
dependence is higher-twist. And secondly, to show that our understanding of higher- 
twist effects as obtained from the 4-pole model, carries over to the BFKL equation. 



10 



Let us start by examining some concrete examples of effective splitting functions. 
Figure 2 illustrates the effective splitting function as a function of x in three situa- 
tions, all with the same value of t, but different sets (a, b and c) of non-perturbative 
parameters, i and tg. Going down in x from x = 1, one sees that initially the three 
splitting functions are almost identical (the inset with the larger scale reveals small 
differences between them). For moderately small x the splitting function actually 
decreases (this phenomenon was recently observed also by [9]), and then starts to 
grow as x~^^. The late onset of the power growth is related to the fact that in cj-space 
the PT pole of the anomalous dimension has a residue of order cSg . 

At a certain point, two of the curves (b and c) change sign (since we use a 
logarithmic scale and plot the absolute value of the splitting function, the change 
of sign appears as a downward cusp) and start to grow with a much larger power 
{uJq). This is the non-perturbative higher- twist component of the splitting function 
discussed in the section 3. As can be seen the exact value of uJq depends on the 
non-perturbative parameters. Furthermore there are situations (curve a) in which 
there is no NP power growth at all, corresponding to the absence of NP zeroes in 
the Green's function. 

Figure 2 is not sufficient to demonstrate that the NP corrections are truly higher 
twist. First we consider the 'PT part' of the splitting function. Figure 3 shows the 
ratio of two effective splitting functions, obtained with different non-perturbative 
parameter sets (a and b), chosen such that there is no component with the large 
NP power growth (which would complicate the interpretation of the ratio). We see 
that the NP parameters affect both the normalisation of the PT splitting function, 
and the exact value of the power growth, since the ratio grows as a power. This 
is as predicted in the 4-pole model, being due to the shift of the PT zeroes of the 
Green's function. We observe that the effect on the power is relatively small (note 
the X scale), and that from a practical point of view it will mostly be the effect on 
the normalisation that will be of interest. The curve at the higher value of t shows 
a significantly decreased dependence on the NP parameters, confirming that their 
effect is higher-twist.^ 

We also need to demonstrate that the NP power component is higher twist. 
Figure 4 shows the splitting function for three different t values, but with the same 
NP parameters. The initial parts of the splitting function now differ since the PT 
scales are different; the NP power growth has (approximately) the same power in the 
three cases, but its normalisation decreases rapidly with increasing t, roughly as e~*, 
confirming that it too is higher twist. 

One of the non-trivial features of the 4-pole model was the (approximate) equality 
between the normalisation of the NP power component, and the quantity u'q — ujq, 



^The attentive reader will have noticed that the modification of the PT exponent changes sign 
at the higher t value — indeed it turns out that the dependence on the NP parameters, while 
decreasing at least as fast as e^*, has a non-trivial t dependence for moderate t values. 
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Figure 2: Examples of effective splitting functions for t = 9.2 (corresponding to Og — 0.12) 
and three different combinations of NP parameters: (a) has i = 1.0, to = 2.0; (b) has 
i = 1.0, to = 3 and (c) has t = 2.0, to = 4.6. The inset shows the same splitting functions 
on an enlarged scale. 

the difference between the exponent of the power growth and the position of the NP 
zero (obtained numerically from the exponent of the NP power growth in the limit 
t -^ oo). Both quantities are shown in figure 5 over a range of t values, illustrating 
very clearly their closeness (as well as their higher-twist nature). A detailed study 
reveals that the relative difference between them is consistent with a term of (ftg), 
as predicted in section 3. We note also that ujq is found to be below the NP exponent 
cup characteristic of the Green's function, consistent with it being due to a zero to 
the left of the leading singularity. 

5. Conclusions 

For some time now there has been some debate as to whether diffusion might destroy 
small-a; factorisation [1] or make it impossible to perturbatively predict the small- a; 
splitting functions [2]. In [4] a model was presented which contained diffusion and 
the correct collinear limits, but also displayed the property of exact factorisation. 
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It failed though to make any statement about (the functional dependence of) the 
magnitude of any higher-twist corrections, leaving open the possibility that they 
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could come to dominate at small x, and thus still destroy factorisation. 

This paper has presented a much more complete study of the problem, both by 
an extension of the model so that it includes leading higher-twist components, and by 
the development of numerical techniques for studying the effective splitting function 
in the full (LL with running coupling) BFKL equation. 

The basic conclusion of these studies is that higher-twist effects, while present 
(and dominant at small x for the effective splitting function) are truly small for 
scaling violations. An explanation of this fact comes from the analysis of the 4-pole 
collinear model, whose Green's function is a sum of two terms, each of which with a 
NP Pomeron part factorised from the t-dependence, the second being uniformly of 
higher twist. 

Although oversimplified, the collinear example is expected to exhibit the mech- 
anism at work in a realistic BFKL equation also. In fact, it can be generalised to 
the case of 2n-poles, for arbitrary n, the Green's function being a sum of n terms of 
higher and higher twist [6]. The basic point remains the fact that the Pomeron sin- 
gularity, although non-perturbative, is stably factorised in front of the t-dependence 
in the leading twist term. 
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A. Green's function of the 4-pole model 

Starting from eq. (3.3) there are various ways of deriving the differential equation of 
the model which differ by the treatment of the higher twist terms. 
The simplest way is to insert the expression 

boot-xil) (A.l) 

into the 7-representation of the (regular) solution 

3^.(t) = e-^*F^(t) = y"^exp [^ - I) t - ^ j\{^') d^' (A.2) 

and to notice that the result vanishes by partial integrations. The corresponding 
differential equation is thus obtained by the replacement 7 — | ^ Sj. By thus using 
the identity 

X{\ + dt) + - = - ., L, , o = -^ (A.3) 



3 df - m + ^ n 



and shifting the t variable to incorporate the constant term 4/3, we obtain the Green's 
function equation 

[uoVt + J\fta,{t)] g{t,to) = [cuA„ +M„as(to)] 9it,to) 

= --K5{t-to) (A.4) 

UJ 

where G^{t,to) = uj~^5{t — to) + e2*-*~*"-'5't^(t, to)Q;s(to) is the gluon Green's function 
discussed in the main text. 

Note first that g{t,to) satisfies a fourth order differential equation in both the 
t and to variables, for any low-t regularisation of as{t), starting from the l/6t ex- 
pression of asymptotic freedom (the regularisation depends on the t parameter of 
the main text). Some care is needed in order to treat the boundary conditions that 
g has to satisfy due to the peculiar distribution occurring in the right-hand side of 
eq. (A.4). This distribution can be taken into account by assuming g{t,to) to be a 
linear combination of left (right) regular solutions of the homogeneous equation, as 
in (3.4), satisfying at t = to the discontinuity requirements 

Ag = Ag" = 0, Ag' = -^^ , A/' = -^ + i^^, (t = to) . (A.5) 
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Given an arbitrary basis jF^(to), J^l(to) of left-regular solutions, the Green's function 
takes the form 

g^it,to) = {:Fl{t):Fl{t,)+T},{t)Tl{t,)))Q{t-t,) 

+ {n{t)n{t^) + n{t)^l{t^))) 0(^0 - t) (A.6) 

and should also be symmetrical under t, to interchange. Because of (A.6), the discon- 
tinuity conditions (A. 5) at t = to can be viewed as a system of four linear equations in 
the four unknowns jF£(to), JF^(to), ^]^(to), ^|;(to), which can be solved by standard 
methods. 

As a consequence, the left-regular solutions JF£, JF| occurring in (3.4) and satis- 
fying (A. 5) are given by the following expressions 

^^^^^ _ detwjAjt), j^ut), nit), nit)) 

^, detiy(^°(t),A(t),^°(t),^l(t)) 

'^^ deiw{n,^hn,n) • ^ ^ 

Here the Ws are the Wronskian matrices of the corresponding functions, where in 
the numerators the column vector of derivatives is replaced, in the proper place, by 
the discontinuity vector (i = 0, 1, 2, 3) 

A(t) = (A,«) = -\ f 0, 4, 0, 7 - 1^^) . (A.8) 

While the Wronskian in the denominator is constant, one can check that the numer- 
ators are indeed solutions of the basic homogeneous equation in (A. 4). Furthermore 
the expressions (A.6) turn out to be independent of the choice of the left-regular ba- 
sis, by the linearity properties of Wronskian matrices. Finally the Green's function 
is determined by symmetry for t < to. 

The explicit determination of J?-]?(to) and jF^(to), and of the corresponding Pomeron 
singularity requires the solution of a matching problem for scattering in a fourth-order 
framework, depending on the regular isat ion of Q;s(t) to the left. The outcome of this 
procedure is the separation of the left-regular solutions into irregular and regular 
ones on the right, as follows 

^°(to) = X^(to) + (Too(cu)^^(to) + ■ ■ ■ , (A.9a) 

nit^) = lUto) + cT,oico)nito) + ■■■ , (A.9b) 

where o"oo(c^), o"io(^)) • • • are scattering coefficients which carry the (non-perturbative) 
Pomeron singularity and X^ ~ i^n)'^, ^r ~ i^K)'^ ^-^^ the irregular solutions on 
the right. The relation of the non-perturbative to perturbative contributions in each 
of the terms of eq. (3.4) is thus very similar to that found in the two-pole model [4]. 
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